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ABSTRACT 

We report on the results of observations of the binary X-ray pulsar LMC X-4 with 
the Rossi X-ray Timing Explorer. Our analysis of the Doppler delays of the 13.5-s X-ray 
pulsations yields the most accurate determination of the LMC X-4 orbital parameters 
available to date. The epoch of orbital phase zero for the 1.4 day orbit is determined 
with an uncertainty of ~ 20 s, and is combined with 5 earlier determinations of the 
epoch of phase zero to obtain the rate of change in the orbital period: P or b/Porb = 
(—9.8 ± 0.7) x 10~ 7 yr~ x . This is the first high significance measurement of the rate of 
change of the orbital period in LMC X-4. 

We present data on one of three strong X-ray flares as well as energy-dependent 
pulse profiles for both non-flaring and flaring time intervals. The pulse profiles during 
the non-flaring time intervals are typically strikingly different from the flare profiles, but 
at other times can be similar. Possible origins of the flaring behavior are discussed. 

We reconsider the orbital decay of LMC X-4 in the context of tidal evolution. We find 
that, while the orbital decay is most likely driven by tidal interactions, the asynchronism 
between the orbit and the rotation of the companion star is probably maintained by the 
evolutionary expansion of the companion star, just as is thought to be the case for 
Cen X-3 and SMC X-l. For LMC X-4, we find that the evidence favors the companion 
star being in a late stage of its evolution on the main sequence. 

The orbital decay timescale of about 1,000,000 years implies that the Roche lobe will 
move catastrophically deep into the atmosphere of the companion within a few times 
10 4 years. This short X-ray lifetime must be similar to the lifetimes of SMC X-l and 
Cen X-3 which have decay timescales of 300,000 and 550,000 years, respectively, and 
may be typical of all massive X-ray binaries in Roche-lobe or near Roche-lobe contact. 

Subject headings: Stars:individual(LMC X-4) — stars:neutron — X-rays:stars 
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1. INTRODUCTION 

The history of the orbital period of a binary star 
system can tell much about the physics of the stellar 
co mponents and their mutual interactions. — In most 



bin aries, the evolution of the orbital period is too slow 

to be detectable, but there are systems in which this 
evolution can be measured. Among X-ray binaries, 
this evolution is apparent in a number of cases. Such 
binaries include Her X-l, which has an intermediate 
mass companion to the compact star, and Cyg X-3, 
which has a companion whose mass is not well deter- 
mined, although it may be a high mass Wolf-Rayet 



star (Deeter et al. 1991, Kitamoto et al. 1995; van 
Kcrkwijk et al. 1996| and references therein). In sys- 



tems such as Cen X-3 and SMC X-l, which comprise 
X-ray pulsars orbiting around high mass companions, 
extraordinary sensitivity to orbital changes and short 
orbital evolutionary timescales combine to allow or- 
bital period changes to be measured with high preci- 
sion. 

Indeed, highly significant detections of orbital pe- 
riod changes have been made in Cen X-3 and SMC X- 
1. Kelley et al. (1983b) determined that the or- 
bital period, P or b, of Cen X-3 is changing at a rate 
Porb/Porb = (-1-8 ± 0.08) x IO7 6 yr-\ while Levine 
et al. (1993) established that P or b/Porb = (—3.36 ± 
0.02) x 10 -6 yr" 1 for SMC X-l. Attempts to measure 
Porb/Porb for LMC X-4 and 4U 1538-52 have resulted 
in only marginal detections or upper limits (see dis- 



cussion in 


Woo et al. 1996 




Safi-Harb, Ogelman, & 


Dennerl 1996], Corbet, Woo, < 


k Nagase 1993). 



Levine et al. (1993) concluded that the rapidly de- 
creasing orbital period in SMC X-l was most likely 
caused by tidal interactions between the neutron star 
and its massive companion. Tidal evolution requires 
asynchronism between the orbital motion and the ro- 
tation of the companion star. The nuclear evolution 
of the companion in its hydrogen shell burning phase 
drives rapid expansion of the star and hence a rapid 
increase in the stellar moment of inertia. Levine et 
al. (1993) argued that this results in a decrease in 
the rate of rotation of the companion star so that 
it would be rotating more 3lowly than required for 



dius is likely to be larger than that of a main-sequence 
star of mass within the range of estimates for Sk 160. 
Tidal interactions most likely also drive the or- 
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bital decay of Cen X-3 (Kelley 
1993). However, Levine et al. (1993) noted that the 
12 ± 2 Rq radius of the companion star in this sys- 
tem is smaller than the 17 ± 1 Rq radius of Sk 160, 
even though both stars have similar masses, and that 
this smaller radius was more likely to be consistent 
with the radius of a ~ 20 Mq star late in its core 
hydrogen-burning phase. Other differences between 
the SMC X-l and Cen X-3 systems were also noted, 
which generally made it difficult to reach firm conclu- 
sions about the exact evolutionary state of the Cen X- 
3 system. 

Safi-Harb, Ogelman, & Dennerl (1996) discussed 
orbital decay of LMC X-4 in the context of conser- 
vative mass transfer and tidal evolution and made a 
comparison between LMC X-4 and Cen X-3. Woo 
et al. (1996) noted that their estimate for the or- 
bital period derivative of LMC X-4, i.e., P or b/Porb = 
(—5.3 ± 2.7) x 10~ 7 yr -1 , could not exclude a fairly 
small value for P or b/Porb- Woo et al. also noted that 
the estimated radius of the companion of LMC X-4, R 
~ 8 Rq, was significantly smaller than that of Sk 160 
even though the masses are not much different, and 
was in the range expected for a ~ 15 Mq star in the 
late stages of its life on the main sequence. Woo et 
al. concluded that the orbital decay rate of LMC X-4 
was consistent with the companion star being in the 
late main-sequence phase. 

LMC X-4 is similar in many respects to Cen X-3 
and SMC X-l. It is a highly luminous X-ray source 
in the Large Magellanic Cloud (Giacconi et al. 1972) 
whose optical counterpart is a 14th magnitude OB 
star dPakull 1976] , |Pcsch, Sandulcak fc Philip 1976| ). 
Pulsations of the X-ray emission with a 13.5 s pe- 
riod demonstrate that the compact object is a neu- 
tron star ( Kelley et al. 198*3^ ) . Optical photometric 
and spectroscopic variations which are periodic and 
X-rays eclipses which recur with the same period of 
1.408 days revealed the binary nature and orbital pe- 
riod of the system (Chevalier & Ilovaisky 1977, Li 



syn chronism with the orbital motion, and would act 

the conditions for orbital decay by tidal torques. The 
Darwin instability was not needed to explain the or- 
bital decay. The inference that hydrogen burning is 
taking place in a shell within Sk 160, the companion 
to SMC X-l, was reinforced by evidence that its ra- 



Rappaport, fc Epstein 197S, White 1978). A 30.5- 



day periodicity in the X-ray intensity, discovered by 
Lang et al. (1981), is most likely due to a precessing 
tilted accretion disk which periodically blocks our line 
of sight to the neutron star. 

Observations of LMC X-4 which we have obtained 
using the Rossi X-ray Timing Explorer (RXTE) are 
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described in Section 2 of this paper. Analyses and 
results pertaining to the timing data, a large flare, 
and pulse profiles are presented in Section 3. The 
results are discussed in Section 4. 

2. OBSERVATIONS 



The R XTE All-Sky Monitor light curve ([Levine et 
ah 1996 ) of LMC X-4 was used to predict the times 
of the LMC X-4 30.5-d cycle high states. A 16-day 
interval, 1998 October 14 to 1998 October 30, was se- 
lected so as to be centered on the time of maximum 
intensity of one of the high states. LMC X-4 was 
then observed with the Proportional Counter Array 
(PCA; Jahoda et al. 1996) and High-Energy X-ray 
Timing Experiment (HEXTE; |Rothschild et al. 1998;) 



on RXTE during 47 intervals within this 16-day time 
interval. All five Proportional Counter Units (PCUs) 
were operating during the great majority of the ob- 
servations, although a few observations were carried 
out with only 3 or 4 operational PCUs. 

The results presented here were obtained from our 
analysis of the PCA data. These data were teleme- 
tered in "Goodxenonl" and "Goodxenon2" modes 
which preserve the 1 microsecond time resolution and 
the inherent energy resolution of the detectors. 

The average count rate attributed to LMC X-4 for 
each observation outside of eclipse is shown in Figure 
1 for each of two energy bands. To make this plot, 
background was subtracted from the observed count 
rates using the faint source version (L7_240) of the 
program "pcabackest." The source intensity varied 
over the 16 day time interval in a manner more or 
less consistent with the general pattern of the 30.5 d 
cycle, i.e, an approximately triangular waveform (see 
Lang et al. 1981 ). The counting rates attributable to 
LMC X-4 reached peak values of ~ 100 and ~ 180 
counts s _1 in the 2-8 and 8-20 keV energy bands, 
respectively, near the expected time of maximum in- 
tensity. We note that the similarity of the variation in 
the two energy bands indicates that the variation can- 
not be attributed solely to photoelectric absorption by 
neutral matter of cosmic element abundances. 

We estimated the X-ray luminosity of LMC X-4 
by integrating a model spectrum consisting of an ex- 
ponentially cut off power law and a Gaussian line at 
~ 6.4 keV derived from a fit to the observed spectrum 
corrected for background. Inclusion of low energy ab- 
sorption in the spectral model generally resulted in a 
negligibly low value for the absorbing column-density 



parameter. The result for the average spectrum ob- 
tained from the observations of 1998 October 21 (Ju- 
lian Date 2451107.85; RXTE observation no. 30085- 
01-15-00) is L ~ 2.3 x 10 38 ergs s" 1 (2 - 25 keV) for 
an assumed distance of 50 kpc. 

3. ANALYSIS AND RESULTS 
3.1. Timing Analysis 

Pulse timing analyses have been performed with 
non-background-subtracted data. The "Goodxenon" 
events were binned into 1/16 s time bins for two en- 
ergy bands, 2-8 keV and 8-20 keV. Only events 
from the front xenon-filled layer of each PCU (i.e., 
LI and Rl) were used for the 2-8 keV band, while 
events from all three xenon-filled layers were used for 
the 8 - 20 keV band. 

For the pulse-timing analysis, the time of each bin 
of data was corrected for the time delay incurred by 
the observation of the pulsations from LMC X-4 at 
the spacecraft position rather than at the barycenter 
of the solar system. In making these corrections we 
used the coordinates of LMC X-4 measured by Bradt, 
Doxsey, & Jernigan (1979) converted to the equinox 
J2000.0 reference frame, viz., a, 72 ooo = 5 h 32 m 49?2, 
<W) = -66° 22' 15". 

Subsequent stages of the pulse timing analysis were 
performed multiple times as part of an iterative pro- 
cess in which the orbital and pulse parameters were 
gradually refined until, as noted below, the changes in 
the parameter values were negligible. Each iteration 
of the procedure began with the further correction 
of the observation times for the orbital motion of the 
pulsar according to a provisional orbital ephemeris for 
LMC X-4. This yields times that would have been 
recorded by an instrument moving with an approxi- 
mately constant velocity relative to the pulsar. The 
corrected times were then used to fold subsets of the 
data according to a provisional pulsar rotation fre- 
quency and frequency derivatives up to second order. 

The count rate data from each of the 47 obser- 
vations were folded to produce pulse profiles. Data 
taken while the source was in eclipse or when a strong 
flare was occurring were not included in the profiles. 
Smearing of these profiles due to inaccuracies in the 
final values of the parameters was less than 0.1 pulse 
periods. 

A pulse template (see Figure 2) was somewhat ar- 
bitrarily constructed by averaging the profiles from 
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the first 12 observations (interval A in Figure 3) which 
yielded high quality phase measurements. The tem- 
plate was then cross correlated with the individual 
profiles, and a precise value for the phase correspond- 
ing to the peak of the cross correlation function was 
obtained via quadratic interpolation of the bin with 
the peak value and the nearest neighbor bins. We ne- 
glected corrections for the zero-delay autocorrelation 
peak that should be present in the cross correlation of 
the template with any of the 12 profiles used to form 
the template, since we believe that such corrections 
are presently unimportant. 

The phases of the cross correlation function peaks 
should be regarded as differences between the ob- 
served pulse phases and those of a model pulsar in 
a model orbit. Phase differences can result from dif- 
ferences between the actual orbit and the model orbit; 
such phase differences can properly be regarded (after 
being multiplied by the pulse period) as pulse arrival 
time differences. Phase differences which result when 
the spin frequency of the actual pulsar differs from 
that of the model pulsar can properly be regarded 
as pulsar rotation phase differences. There can also 
be other causes of phase differences, such as changes 
in the intrinsic beam pattern of the pulsar or simply 
statistical fluctuations in the observed count rates. 

We performed the pulse timing analysis on the data 
from each of the two energy bands separately. The 
cross correlation phase differences with respect to a 
best-fit model described below are shown in Figure 3. 
This figure clearly shows that phase differences which 
are small in magnitude, i.e., < 0.1 cycles, and rela- 
tively constant from one observation to the next were 
obtained for a majority but not for all of the obser- 
vations. The large scatter of the phase differences for 
the first five observations, which were obtained when 
the overall source intensity was low, cannot be en- 
tirely attributed to lack of sufficient signal. Rather, 
the pulse profile was significantly different from the 
template profile during three of the five observations, 
while the pulsations were very weak during two of the 
observations. 

The pulse phases obtained from the last ~ 3 days 
of the observations are not consistent between the re- 
sults for the two energy bands and also show both sig- 
nificant changes (relative to our best fit model) as well 
as sudden jumps. These problems are mostly caused 
by substantial changes in the pulse profile which are 
apparent by direct inspection and by examination of 
the relative strengths of the harmonics of the pulse 



frequency in Fourier transforms. The pulse profile 
changes are discussed further in Section 3.3. 

We fit the results for the 28 observation intervals 
in which the phase differences appeared to be well- 
behaved (solid black circles in Figure 3) with a general 
circular orbit model to obtain corrections to the pro- 
visional orbital and pulse period ephemerides. The 
model included parameters for the epoch of orbital 
phase zero, T^/ 2 , and projected orbital radius, a x sini, 
as well as parameters for the pulse epoch, frequency, 
and two frequency derivatives. Corrections to model 
parameters were obtained using a linearized (with 
respect to a provisional model) unweighted least- 
squares fitting procedure. The root mean square de- 
viation of the observed arrival times from the model 
arrival times was used as an estimate of the obser- 
vational uncertainty in each arrival time. The latter 
uncertainties were then propagated through the linear 
least-squares procedure to estimate la errors in the 
orbital and spin frequency parameters. The folding 
and fitting procedures were iterated using corrected 
parameters until the corrections became small in com- 
parison with their estimated uncertainties. 

The parameters for the best fit circular orbit model 
are given in Table 1. We derive a value for a^sini 
of 26.34 ± 0.02(lcr) lt-s which is consistent with the 
value of 26.31 ± 0.07 lt-s derived from Ginga obser- 
vations by Levine et al. (1991), with the value of 
26.27 ± 0.04 lt-s derived from an earlier Ginga ob- 
servation by Woo et al. (1996), and with slightly 
less precise values derived from ROSAT data by Safi- 
Harb, Ogclman, & Dennerl (1996) and Woo et al. 
(1996). The weighted average of the fitted values of 
the epoch of orbital phase zero for the two energy 
bands is T n/2 = JD 2, 451, 1 11. 86579 ± 0.00010 (TT). 
The estimated error in T n / 2 is calculated assuming 
statistical independence of the errors in the results 
from each energy band. These errors are, in turn, 
estimated on the basis of the scatter of the observed 
arrival times relative to the model arrival times as 
described above. However, the difference between the 
values of T n / 2 obtained for the two energy bands indi- 
cates that the magnitude of the error of the weighted 
average may be underestimated by a factor as large 
as ~ 2, and we therefore use an error estimate of 
0.00020 d (= 17 s)in our determination of the orbital 
decay (see below). In any case, this result is signifi- 
cantly more accurate than previous determinations of 
the epoch of phase zero of LMC X-4. 

We have also fit the pulse timing results from the 
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28 observations with a model that allows a small ec- 



centricity e by including parameters that are propor- 



tional to e sin Q and e cos f2 as coefficients of cos 2ui or bt 
and sin2w or ;,t, respectively. Here, il is the longitude 
of periastron and Lu rb is the orbital angular frequency. 



No significant eccentricity was found, so that we re- 



port an upper limit e < 0.003 (2a). 

The resultant Doppler delay data, with the delays 
due to the best-fit circular orbit added back in, are 
plotted in Figure 4. 

Our determination of the epoch of orbital phase 
zero is shown together with previously determined 
epochs of orbital phase zero for LMC X-4 in Fig- 
ure 5. The solid curve is the best-fit quadratic 
function of time (see Table 1). The coefficient of 
the quadratic term of the fitting function implies 



Porbj Pq- 



(-9.8 ± 0.7) x 10" 7 yr" 1 . This repre- 



sents the first highly significant evidence for orbital 
decay in LMC X-4. It is consistent with the esti- 
mate of Woo et al. (1996) who obtained P or b/Porb = 
(-5.3 ±2.7) x 10~ 7 yr" 1 . The coefficients for the best 
fit quadratic function that gives the predicted times 
of orbital phase zero in LMC X-4 may be found in 
Table 1. 

We obtain a value for the pulse period (Table 
1) which is barely shorter than that measured in 



19 91 with ROSAT (13.50292 ± 0.00002 s; 3afi-Harb 
Og fclman, fc Denncrl 199^ , |Woo ct al. 1996| ). The a7 



erage value of the pulse period derivative over the 
intervening ~ 7 years is P pu ise/ Ppuise = (-4.9 ± 



0.2) x 10" 



yr~ 



which is much smaller in mag- 



nitude than the derivative which we estimate from 
our timing analysis of the RXTE observations, i.e., 
Ppuise/ Ppuise = (-2.14 ±0.02) x 10~ 3 yr^ 1 . A graph 
of the history of the pulse period in Woo et al. (1996) 
shows that episodes of both spin up and spin down 
have occurred during the last 20 years. Since the long- 
term average spin up (or spin down) rate of LMC X-4 
is much smaller in magnitude than the rate typically 
measured during an observation of a few weeks or less 
in duration, the pulsar may be near its equilibrium ro- 
tation period (see also Naranan et al. 1985] , Woo et 
al. [1996]). 



3.2. Large Flares 

It is believed that the X-ray luminosity of LMC 
X-4 is usually fairly constant at a value somewhat 
higher than the Eddington luminosity for a 1.4 Mq 
neutron star, even though the X-ray flux observed at 



Earth is modulated with a 30.5-day period (Lang et 
al. 1981). It has been known for a long time, how- 



ever, that LMC X-4 exhibits large X-ray flares, which 
typically last for ~ 1000 sec and which occur non- 



pcriodically approximately every few days (Kelley et 



al. 1983a, Levine et al. 1991 and references therein). 



During the flares, the 2-25 keV X-ray luminosity in- 
creases by factors of ~ 2 to 5, and the X-ray spectrum 
softens considerably. 

Flares were apparent during the RXTE observa- 
tions of 1998 October 14 (Julian Date ~ 2451100.9, 
RXTE observation no. 30085-02-01-00), October 26 
(~ 2451112.9, 30085-01-27-02), and October 28 (~ 
2451114.1, 30085-01-29-00). A plot of the count rate 
during a portion of the largest of these flares (that 
of 1998 October 28) is shown in Figure 6. The ob- 
servations were interrupted just after a peak in the 
flare by an Earth occultation of the source as seen 
from the moving RXTE spacecraft. The top panels 
of Figure 6 show the count rate of LMC X-4 in 2- 
sec time bins for the 8-20 keV and 2-8 keV energy 
bands. Three temporal structures are evident in the 
raw counting rate: the 13.5-sec pulsations, a modula- 
tion with a characteristic timescale of ~ 150 sec, and 
an overall outburst time of > 400 sec. The bottom 
panel of Figure 6 shows the ratio of the count rate 
in the 2-8 keV to that in the 8-20 keV energy chan- 
nel, i.e., a "softness" ratio. As was found in previous 
observations of these flares, they are spectrally very 
soft. Note also that the pulsations essentially disap- 
pear in the softness plot, indicating that the softness 
is a monotonically increasing function of the source 
luminosity and is not a function of pulse phase. 

The peaks of the X-ray counting rates during this 
particular flare are up, compared with the quiescent 
level, by factors of 4 and 7 in the 8-20 keV and 2-8 keV 
bands, respectively. After properly integrating over 
the 2-25 keV X-ray band and averaging over at least 
one pulse cycle, we find that the X-ray luminosity 
increased by a factor of ~ 3, which yields an absolute 
value for the peak luminosity of the flare of ~ 6 x 10 38 
ergs s^ 1 . 

We discuss physical implications of the large X-ray 
flares in section 4.2. 

3.3. Pulse Profiles 

We used the best-fit orbital and pulse frequency 
parameters to produce three sets of pulse profiles for 
each of 10 energy intervals (Figures 7-9). The profile 
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sets are averages for (1) the same 12 observations used 
to form the template pulse profile (interval A of Fig. 
3), (2) the portion of the observation of 1998 October 
28 during which a strong flare was observed (interval 
B) , and (3) three of the later observations (interval C) 
in which the 8-20 keV pulse phase was quite different 
from that anticipated from the 28 observations used 
for the pulse timing fits. 

The profiles in Figure 7 can be compared with the 
non-flare profiles obtained from Ginga observations in 
1988 (shown in |Woo et al. 1996Q and 1989 (shown in 
Levine et al. 1991). As a rough approximation, the 



profiles can be thought of as consisting of an under- 
lying sinusoid with phase and strength that may vary 
with photon energy and a prominent pair of narrow 
absorption dips. In both the Ginga and RXTE sets of 
non-flare profiles the narrow dips (at phases ~ 0.58 
and 0.82 in Figure 7) are separated by Atfi ~ 0.25 
and the first dip is more prominent, i.e., deeper, than 
the second. In general, the dips exhibit structure on 
time scales down to the time resolution of the folded 
profiles, e.g., ~ 0.27 s in Figure 7. In the Ginga pro- 
files the underlying sinusoidal profile appears to be 
stronger at low energies and shifted in phase (rela- 
tive to the dips) in comparison with the underlying 
sinusoid in the present RXTE profiles. 

Average pulse profiles during the time of the strong 
flare are shown in Figure 8. The epoch of pulse phase 
zero is the same as for the profiles shown in Figures 7 
and 9. In general, the pulse profiles during the flare 
are simpler than those outside of the flare and are es- 
sentially sinusoidal at all energies studied. Moreover, 
during the flare the pulses are in phase at all energies, 
in contrast with the profiles outside of the flares. We 
also note that the flare profiles are very similar to the 
flare pulse profiles obtained from Ginga observations 
in 1989 ( 



Levine et al. 1991) 



The set of pulse profiles produced from data taken 
late in the RXTE observation of LMC X-4, and during 
an interval when there was no obvious flaring activity, 
is shown in Figure 9. These profiles are very similar 
to those obtained during the flare (see Figure 8) in 
that they are rather sinusoidal in shape, and in phase 
at all energies. However, the fractional amplitude of 
the modulation is much lower in the profiles in Fig. 9 
than in the flare profiles. This is a dramatic example 
of how the pulse profiles can change long (~ a day) 
after the occurrence of a major flare, although it is 
possible that strong flares may have occurred closer 
to these later observations during the large gaps in 



observational coverage. It is changes in pulse shape 
such as these that account for the occasional erratic 
behavior of the pulse phase, as exhibited in Figure 3. 

4. DISCUSSION 

4.1. Orbital Period Changes 

Levine et al. (1993) discuss orbital period changes 
in the context of the measured orbital period deriva- 
tive of SMC X-l. They derive an estimate for the 
orbital period derivative for the situation where ex- 
pansion of the companion star drives the orbital de- 
cay. Their discussion is pertinent when the logarith- 
mic expansion rate of the companion is roughly con- 
stant for at least the time necessary to establish a 
pseudo-equilibrium in the difference between the ro- 
tation frequency of the companion star and the orbital 
frequency. In this c ase the orbital peri od derivative is 
given by (eqn. 6 of Levine et al. 1992) 



orb 
°orb 



lu c dln(I)/dt 



(1) 



where lu c and lok are, respectively, the angular fre- 
quencies of the companion's rotation and the orbital 
motion, / is the moment of inertia of the companion, 
fi is the system reduced mass, and a is the separation 
of the centers of the two stars. 

Monte Carlo analysis of the binary system param- 
eters of LMC X-4 yields an estimated mass of the 
companion of 12 to 18 Mq (see Joss & Rappaport 



1984 |Nagase 1989] , |Levine et al. 1991| ). Stellar evolu 
tion calculations for a 15 Mq star (Ph. Podsiadlowski 
1991, private communication) show that the rate of 
increase of the companion star's moment of inertia 
can be as large as dln(I)/dt ~ 3 x 10~ 7 yr _1 near 
the end of the main-sequence phase when the sup- 
ply of hydrogen in the star's core is nearly exhausted. 
The Monte Carlo analysis of the system parameters 
also indicates that the factor (/xa 2 /3/— 1) most likely 
has a value in the range of 0.3 to 1.1. This would 
then yield an orbital period derivative P or b/Porb ~ 
— (3 — 10) x 10~ 7 yr _1 which is compatible with the 
measured value Porb/Porb = {— 9. 8 ±0.7) x 10~ 7 yr _1 . 

The radius of the companion star is estimated from 
t he Monte Carlo analysis to be in the r a nge 6-9 R^ 



Joss fc Rappaport 1984, Nagase 1989, Levine et al 



1991) which is compatible with the radius near the 
end of the main-sequence phase derived by Podsiad- 
lowski. 



G 



One problem is that the optically-observed lumi- 
nosity and temperature place the companion in a re- 
gion of the H R diagram which indicates a mass of 
~ 25 Mq, which is above the range indicated by the 
Monte Carlo calculations. We note that the calcula- 
tions of Podsiadlowski did not incorporate several ef- 
fects which could be important in the determination 
of the companion mass, e.g., mass loss via a wind, ef- 
fects of the evolution of the companion in a binary sys- 
tem with the likelihood that substantial mass transfer 
has occurred, the low metallicity of stars in the LMC, 
and effects of X-ray illumination. 

Our measured orbital period derivative also firmly 
establishes that the orbit of LMC X-4 is decaying, as 
are the orbits of SMC X-l and Cen X-3, on a timescale 
of about half a million years. Since the photospheric 
gas-pressure scale height of the massive companion 
star is expected to typically be about 0.1i?Q, then 
as the orbit decays the Roche lobe will descend into 
the companion's atmosphere by an additional scale 
height every ~10,000 yr. Thus, the mass transfer will 
run away to super-Eddington rates over a relatively 
short timescale. This process will presumably end 
with the neutron star spiraling into the envelope of 
the companion, thereby terminating the high-mass X- 
ray binary phase of the evolution. See Levine et al. 
(1993) for an extensive discussion of this scenario in 
the context of SMC X-l. 

4.2. Large X-Ray Flares 

The three large flares seen in the present observa- 
tions are similar in all respects to the flares of LMC X- 



4 observed in detail on earlier occasions (Levine et 



al 1991 and references therein). During such flares 
the X-ray luminosity (2-25 keV) goes from its typical 
nearly steady value of ~ 2 x 10 38 ergs s" 1 up to ~ 10 39 
ergs s _1 , and the modulation factor of the pulsations 
increases substantially. As the luminosity (computed 
using the flux averaged over a pulse cycle) increases 
by factors of 2 to 5, the X-ray spectrum softens con- 



sid erably. While the X-ray spectrum during neither 
the flaring nor quiescent modes is well fit by a black- 
body, we can say that the characteristic temperature 
associated with the emission decreases by about a fac- 
tor of 2 during the flaring state. This was shown for 
the instance of the flare observed on 1998 October 
28 by fitting the spectrum with the exponential form 
e ~E/kT _ rp ne -(.Q-j-gj ener gy released in a flare is typi- 
cally ~ 10 41 ergs. Finally, we point out that no other 
known accretion-powered X-ray pulsar exhibits such 



flares. 

From the perspective of a simple thermalization 
scenario, a factor of 2 to 5 increase in luminosity, 
coupled with a decrease in temperature by a factor 
of ~ 2, leads us to infer that the emitting area in- 
creases by a factor of about 50! For a nominal polar 
cap surface area of ~ 3 km 2 , which might correspond 
to the accreting area during the quiescent state, the 
emitting area during the flare would have to increase 
to more than ~ 100 km 2 . This is equivalent to the 
surface area of a cylindrical accretion column 25 km 
tall with a diameter of 2 km, or of a cap that extends 
~ 40° outward from the pole. 

The obvious question that arises is "what is the 
physical origin of the large flaring events?" We dis- 
cuss briefly three possible origins for the flares: (1) 
an instability which dumps ~ 10 21 grams of matter 
stored in the accretion disk or magnetosphere, with 
the attendant release of gravitational potential en- 
ergy; (2) the release of nuclear energy due to the fu- 
sion of carbon and similar elements below the surface 
of the neutron star; and (3) an event involving the 
release of magnetic energy. 

It is commonly thought that variations in accretion 
rate are responsible for the ubiquitous variability and 
the flaring seen in a large fraction of X-ray binaries. 
Accretion events must therefore be considered as can- 
didate sources of energy for the LMC X-4 flares. It 
is of interest to note that the energy released by the 
Type II bursts that are attributed to accretion in- 
stabilities in GRO J 1744-28 can be comparable with 
that released by the flares in LMC X-4; furthermore, 
the GRO J1744-28 bursts may exceed the Edding- 
ton luminosity for a canonical neutron star by factors 
up to - 100 ( |Jahoda et al. 1998j ). On the other 
hand, the flaring in LMC X-4 is quite distinct in de- 
tail from the variability and flaring in other X-ray 
binaries, including the Type II X-ray bursts in the 
Rapid Burster (MXB 1730-335) and GRO J1744-28 
(pucrriero et al. 199S| , [Nishiuchi et al. 1999| , [Woods et 



al. 1999 and references therein). Moreover, in a high- 



strength magnetic field pulsar such as LMC X-4, the 
inferred increase in surface area during the flares does 
not seem to have a ready physical interpretation. For 
example, emission from the sides of an accretion col- 
umn might require a height of more than a neutron 
star radius. In order to confidently retain the picture 
of the flares as accretion events, one would like to be 
able to model accretion onto the polar caps of a neu- 
tron star at the Eddington limit, increase the rate by 
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a factor of ~ 5 for ~ 500 sec, and compute the ex- 
pected spectrum at least well enough to demonstrate 
the observed softening. To our knowledge, this has 
not been done. 

An energy release of ~ 10 41 ergs, within the con- 



text of a nuclear burning scenario (Brown & Bildstcn 
19 d|) would require burning ~ 10 23 g of carbon and 



released. If this is true, then the field would need to 
be regenerated, since observations suggest that flares 
have occurred every few days over the past 20 years, 
at least. To our knowledge, the detailed physics of 
such a scenario has not yet been explored. 

A number of similarities between LMC X-4 and 



soft gamma repeaters (SGRs, Kouveliotou 1999 



similar elements, with a nuclear burning efficiency of 



10 



(Hydrogen and helium are burned almost 



Hur 



ley 2000) suggest a possible connection within the 



context of the magnetar model (Thompson & Dun- 



as quickly as they are accreted very near the surface 



of the neutron star.) This scenario has the advan- 
tage that if the nuclear burning takes place over a 
region of the neutron star surface that is larger than 
the footprint of the accretion flow, the flare spectrum 
would naturally be expected to be softer than the 
quiescent spectrum. At accretion rates of 10 18 g s _1 , 
the mean time between outbursts would be ~ 10 5 s 



can 1995). First, even though the large flares in LMC 
X-4 last much longer than the bursts from SGRs, the 
energies of the flares and bursts are roughly compa- 
rable; the flares in LMC X-4 have fluences of ~ 10 41 
ergs and SGR bursts have fluences often in the range 
10 38 to 10 41 ergs QGogiis. et al. 200q). Second, if the 



neutron star in LMC X-4 is rotating at about the 
equilibrium rate estimated in the standard accretion 



(~ 


1 day) , which is roughly typical of flare recurrence 


torque model ( 


Lamb, Pethick, & Pines 1973 




Rap- 


tim 


es. The large and fairly constant ~ 50% mod- 


paport & Joss 


1977 




Ghosh & Lamb 1979 


. Bildsten 


ula 


;ion factor of the pulses during the flares, com- 


ct al. 1997), this could imply a surface dipole field 



pared with only ~ 20% away from the flares, might 
be explained by subsurface nuclear burning that takes 
place at only one of the two magnetic poles during 
a particular event, so that substantial modulation 
would be observed as the polar cap rotates in and 
out of view. However, the subflare structure observed 
on a timescale of ~ 150 s lacks a natural explanation 
in the context of a deep nuclear burning scenario. Fi- 
nally, in this regard we note that Brown & Bildsten 
(1998) found difficulties explaining the large flares in 
LMC X-4 with carbon burning. The principal diffi- 
culty is the fact that accretion of H-rich material is 
thought to produce much heavier "ashes" than car- 
bon, in which case not much carbon would be ex- 
pected to accumulate and be available to power the 
flares. 

Finally, we consider the possibility that magnetic 
energy powers the flares. Magnetic energy might be 
suddenly released, e.g., if the accreting material at the 
polar caps were somehow able to compress the mag- 
netic field until an instability leads to a reconnection 
of the field lines. The energy of ~ 10 41 ergs in a flare 
would require that the product of magnetic field and 
volume be 

Bf 3 V km3 ~ 25 (2) 

where -B13 is the surface magnetic field strength in 
units of 10 13 G, and V^ m 3 is the volume in units of 
cubic kilometers. For B\% = 1, the magnetic energy 
contained within ~ 1% of the neutron star volume or, 
alternatively, of its magnetosphere would need to be 



strength of ~ 3 x 10 13 G. This is stronger than the 
canonical 10 12 G field strength thought to be typi- 
cal of accretion-powered pulsars, but is somewhat be- 
low the field strengths thought to characterize mag- 
netars. However, we note that Makishima et al. 
(1999) estimate a surface field strength for LMC X- 
4 of ~ 2 x 10 12 G based on the shape of its high 
energy spectrum. Third, both LMC X-4 and the 
persistent pulsing sources associated with SGRs have 
notable variations in pulse profile shapes (C. Kou- 
veliotou 2000, talk given at Rossi 2000, Greenbelt, 
MD). Fourth, two of the five known SGRs are nearly 
coincident in direction with clusters of massive stars 



flFuchs et al. 1999] , |Vrba et al. 2000] ), which suggests 
that SGRs might be associated with massive stars. 
Perhaps the facts that the LMC X-4 flares are longer 
in duration and have softer spectra than the bursts 
of soft-gamma repeaters is related to a lower mag- 
netic field strength and the presence of accretion in 
LMC X-4. 

Thus, at this point, we are uncertain as to whether 
the origin of the energy in the flares is gravitational, 
nuclear, or magnetic. However, LMC X-4 is unique in 
its flaring behavior, and the resolution of the origin 
of the flares may well provide a key to understanding 
some fundamental property of the neutron star. 



The authors gratefully acknowledge useful conver- 
sations with Lars Bildsten, Andrei Gruzinov, Ed Mor- 
gan, Mike Muno, and Henk Spruit. We also thank 



8 



Mike Muno for assistance with the spectral analysis. 
This work was supported in part by NASA Grant 
NAG5-8244 and by NASA Contract NAS5-30612. 

REFERENCES 

Bradt, H. V., Doxsey, R. E., & Jernigan, J. G. 1979, 
in X-ray Astronomy, ed. W. A. Baity & L. E. Pe- 
terson, p. 3, (Oxford:Pergamon) 

Bildsten, L., et al. 1997, ApJS, 113, 367 

Brown, E. F., & Bildsten, L. 1998, ApJ, 496, 915 

Chevalier, C, & Ilovaisky, S. A. 1977, A&A, 59, L9 

Corbet, R. H. D., Woo, J. W., & Nagase, F. 1993, 
A&A, 276, 52 

Deeter, J. E., Boynton, P. E., Miyamoto, S., Kita- 
moto, S., Nagase, F., & Kawai, N. 1991, ApJ, 383, 
324 

Dennerl, K. 1991, Thesis - Ludwig-Maximilians Univ. 
Max-Planck-Inst. fur Physik und Astrophysik, 
Garching (Germany) 

Fuchs, Y., Mirabel, F., Chaty, S., Claret, A., Ce- 
sarsky, C. J., & Cesarsky, D. A. 1999, A&A, 350, 
891 

Ghosh, P., & Lamb, F. K. 1979, ApJ, 234, 296 

Giacconi, R., Murray, S., Gursky, H., Kellogg, E., 
Schreier, E., & Tananbaum, H. 1972, ApJ, 178, 
281 

G6gii§, E., Woods, P. M., Kouveliotou, C, van 
Paradijs, J., Briggs, M. S., Duncan, R. C, 



& Thompson, C. 2000, ApJ, in press (astro- 
])h/0002181|) 



Guerriero, R., et al. 1999, MNRAS, 307, 179 

Hurley, K. 2000, in Proceedings of the 5th Huntsville 
GRB Symposium (Melville: AIP Press), in press 
(|astro-ph/9912061|) 



Kelley, R. L., Jernigan, J. G., Levine, A., Petro, L. 
D., & Rappaport, S. 1983a, ApJ, 264, 568 

Kelley, R. L., Rappaport, S., Clark, G. W., & Petro, 
L. D. 1983b, ApJ, 268, 790 

Kelley, R. L. 1986, in The Evolution of Galactic X-ray 
Sources, ed. J. Truemper et al. (Dordrecht:Reidel), 
75 



Kitamoto, S., et al. 1995, PASJ, 47, 233 

Kouveliotou, C. 1999, Proc. Natl. Acad. Sci., 96, 5351 

Jahoda, K., Swank, J. H., Giles, A. B., Stark, M. J., 
Strohmayer, T., Zhang, W., & Morgan, E. 1996, 
Proc. SPIE, 2808, 59 

Jahoda, K., Stark, M. J., Strohmayer, T. E., Zhang, 
W., Morgan, E. H., & Fox, D. 1998, Nucl. Physics 
B (Proc. Suppl.), 69/1-3, 210 

Joss, P. C, & Rappaport, S. 1984, ARA&A, 22, 537 

Lamb, F. K., Pethick, C. J., & Pines, D. 1973, ApJ, 
184, 271 

Lang, F. L., et al. 1981, ApJ, 246, L21 

Levine, A. M., Bradt, H., Cui, W., Jernigan, J. G., 
Morgan, E. H., Remillard, R., Shirey, R. E., & 
Smith, D. A. 1996, ApJ, 469, L33 

Levine, A., Rappaport, S., Deeter, J. E., Boynton, P. 
E., & Nagase, F. 1993, ApJ, 410, 328 

Levine, A., Rappaport, S., Putney, A., Corbet, R., & 
Nagase, F. 1991, ApJ, 381, 101 

Li, F., Rappaport, S., & Epstein, A. 1978, Nature, 
271, 37 

Makishima, K., Mihara, T., Nagase, F., & Tanaka, Y. 
1999, ApJ, 525, 978 

Nagase, F. 1989, PASJ, 41, 1 

Naranan, S., et al. 1985, ApJ, 290, 487 

Nishiuchi, M., et al. 1999, ApJ, 517, 436 

Pakull, M. W. 1976, IAU Circ, No. 3017 

Rappaport, S., & Joss, P. C. 1977, Nature, 266, 683 

Rothschild, R. E., et al. 1998, ApJ, 496, 538 

Safi-Harb, S., Ogelman, H., & Dennerl, K. 1996, ApJ, 
456, L37 

Pesch, P., Sanduleak, N., & Philip, A. G. D. 1976, 
IAU Circ, No. 3023 

Thompson, C, & Duncan, R. C. 1995, MNRAS, 275, 
255 

van Kerkwijk, M. H., Geballe, T. R., King, D. L., van 
der Klis, M., & van Paradijs, J. 1996, A&A, 314, 
521 



9 



Table 1 

Orbital and Pulse Parameters for LMC X-4 



Parameter 


Units 


Value a (2-8 keV) 


Value a (8-20 keV) 


a x sini b 


lt-s 


26.333 ±0.019 


26.370 ±0.031 


e 




< 0.003 (2a) 




rp C 

J-tt/2 


MJD d TT 


51110.86571 ±0.00012 


51110.86600± 0.00020 


^ 'pulse ' 


Hz 


0.074060687 ± 0.000000010 


0.074060697 ± 0.000000016 


V -pulse. 


Hz s" 1 


(5.01 ±0.04) x 10~ 12 


(5.11 ±0.07) x 10~ 12 


V pulse 


Hz s~ 2 


(-1.4 ±0.7) x 10~ 18 


(-1.2 ±1.2) x 10~ 18 


a s 


MJD TT 


48, 137.7500 ±0.0004 




ai s 


days 


1.40839776 ±0.00000026 




a 2 g 


days 


(-2.65 ±0.19) x 10~ 9 




Porb J Porb 


yi-- 1 


(-9.8 ±0.7) x 10~ 7 





a The errors (la) have been obtained from a least-squares analysis in which the la- 
errors of the observed pulse arrival times were estimated by the root mean square 
deviation of the measurements from the best fit model. This may underestimate the 
effects of certain types of systematic errors. 

h a x sin? weighted average = 26.343 ± 0.016 lt-s 

c T w/2 average = MJD 51110.86579 ± 0.00010 (TT) (but see text) 

d Modificd Julian Date = Julian Date - 2,400,000.5 

Tor the epoch MJD 51235.0 

{ P P ulse = 13.502440 ± 0.000002 s (2-8 keV) or 13.502439 ± 0.000003 s (8-20 keV) 
s The predicted times of orbital phase zero are given by: = a + a\N + GI2-/V 2 . 
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Fig. 1. — Average background-subtracted count rates 
in 2 energy bands for each of the 47 observations of 
LMC X-4 obtained in 1998 October with the PCA on 
RXTE . The observations in which strong flares oc- 
curred are denoted by "*" ; however, only the data 
taken during non-flaring intervals were used in com- 
puting the average rates shown for these observations. 
Time intervals during which LMC X-4 was in eclipse 
were also excluded in the calculation of the average 
rates. In the few instances when fewer than 5 PCUs 
were in operation, the rates have been normalized to 
a 5 PCU basis. Estimated errors in the average rates 
due to counting statistics are negligible. Julian Date 
2,451,100.5 corresponds to 1998 October 14 h TT. 



This 2-column preprint was prepared with the AAS IATjrjX 
macros v4.0. 
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Fig. 2. — Average pulse profiles for 2 energy bands 
for 12 observations of LMC X-4 (see interval A in 
Fig. 3). These profiles were used as the templates for 
our pulse timing analysis. The la uncertainties from 
counting statistics are approximately 0.5 counts s _1 . 
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Fig. 3. — Results of our pulse timing analysis for 2 
energy bands. Each symbol gives the phase of the 
pulse profile for one of the 47 observations of LMC 
X-4 relative to the phase of a model pulsar. The solid 
circles indicate measurements that we use in the least- 
squares fit to obtain orbital and spin parameters. The 
open circles are measurements that are not used in 
the fit because the timing results here indicate signif- 
icant inconsistencies from the results expected from a 
pulsar with an invariant pulse profile. The observa- 
tions in which strong flares occurred are denoted by 
"*" ; however, the data taken during the flares were 
not included in the pulse profiles made for the timing 
analysis. Time interval A was used for the construc- 
tion of the template pulse profiles (Fig. 2). Time 
intervals A, B, and C were used to make the energy- 
dependent profiles shown in Figures 7-9, respectively. 
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Fig. 4. — Inferred pulse arrival time delays for 28 
observations of LMC X-4 (see Fig, 3 and text). The 
solid circles and sine wave fit show the delays with 
respect to a model pulsar at the center of mass of the 
binary system (scale on left), while the crosses show 
the measured time delays relative to a pulsar in the 
best fit orbit (scale on right). 



Fig. 5. — Epochs of orbital phase zero of LMC X-4 
are shown relative to the epochs expected according 
to the best- fit constant orbital period of 1.40839457 
days. The measurements are based on EXOSAT, 
Ginga, and ROSAT data (Dennerl 1991; Levine et 
al. 1991; Woo et al. 1996) and the present work. The 
curve shows the best-fit quadratic function. The error 
bars shown indicate ±lcr uncertainties. 
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200 300 400 

Time (s) 

Fig. 6. — (top and middle) Background-subtracted 
count rates for 2 energy bands during a portion of 
the observation of 1998 October 28 (RXTE observa- 
tion no. 30085-01-29-00) when a strong flare was in 
progress. See Section 3.2 for a detailed description of 
the flare. The observation was terminated by Earth 
occultation of the source at about 430 s while the flare 
was in progress. The 13.5 s pulsations arc highly ev- 
ident in both energy bands, (bottom) The softness 
ratio defined by the number of counts in the 2-8 keV 
band divided by the number in the 8-20 keV band. 
The rates and the softness ratio are averaged in 2 s 
time bins. The 13.5 s pulsations are not evident in 
the softness ratio plot. 
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Fig. 7. — Average background-subtracted pulse pro- 
files for 10 energy bands for the 12 observations in 
time interval A (see Fig. 3). Each profile is plotted 
twice for clarity. A typical ±1ct error and the energy 
range in keV are shown for each profile. 
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Fig. 8. — Average background-subtracted pulse pro- 
files for 10 energy bands for the observation of the 
flare which occurred on 1998 October 28 (during in- 
terval B of Fig. 3). The profiles were constructed 
using data from the entire time interval plotted in 
Fig. 6. Each profile is plotted twice for clarity. A 
typical ±1(7 error and the energy range in keV are 
shown for each profile. 
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Fig. 9. — Average background-subtracted pulse pro- 
files for 10 energy bands for the 3 observations in time 
interval C (see Fig. 3). Each profile is plotted twice 
for clarity. A typical ±1<7 error and the energy range 
in keV are shown for each profile. 
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